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ABSTRACT 


We show how, as a result of the strong heating produced at chromospheric 
levels during a solar flai . burst, the local gas pressure can transiently 
attain very large values in certain regions. The effectiveness of the sur- 
rounding magnetic field at confining this hi^ pressure plasma is therefore 
reduced and the flaring loop becomes free to expand laterally. In so doing it 
may drive magnetic field lines into neighboring, non-flaring, loops in the 
same active region, causing magnetic reconnection to take place and triggering 
another flare burst. The features of this interacting loop model are found to 
be in good agreement with the energetics and time structure of flare- 
associated solar hard X-ray bursts. 

Subject headings ; hydroma gne tics - plasmas - Sun ; flares - Sun : X-rays 
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!• INTRODUCTION 

Solar hard X-ray bursts (photon anargy e X 10 kaV) fraquantly axhibit da- 
tallad time structura on tiniascalas of a few saconds (saa, a*g.^ Hoyng, Brown, 
and van Baek 1976; Dannis, Frost, and Orwig 1981). This has lad to tha con- 
capt of tha "Elementary Flare Burst" (EFB - van Beak, de Feiter, and da Jagar 
1974; da Jager and da Jonga 1978), whereby it is supposed that hard X-ray 
bursts with a ccMnplicatad time structure involving several "spikes" are in 
fact c<Mi^8ed of a number of discrete EFB's, each with a simple "rise-fall" 
time structure. 

de Jager and de Jonge (1978) were the first to point out that there ap- 
pears to be a significant difference (as measured by their FWHMs) amongst 
EFB's in a particular hard X-ray event, suggesting that different EFB's cor- 
respond to the successive activation of a number of different sources, rather 
than the continual reactivation of a single source. Karpen, Crannell, and 
Frost (1979) also studied the structure of hard X-ray bursts exhibiting clear 
EFB structure (in their nomenclature, "multiply impulsive" events), together 
with microwave data on the same bursts, and showed that the source parameters 
(density, temperature, magnetic field, etc.) for each EFB were markedly dif- 
ferent. This similarly led them to the conclusion that different EFB's origi- 
nate in different regions of the flare. Finally, Kane, Pick, and Raoult 
(1980) studied the spatial structure of type III radio bursts observed to be 
synchronous with hard X-ray EFB's and deduced that each of the hard X-ray 
bursts in the event in question (on 2 September 1978) originated in one of two 
spatially distinct regions, separated by some 7 minutes of arc, or 3 x 10^® 


cm . 
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Th^re thus appears to be a significant anount of evidence for different EFB’s 
corresponding to excitation of different parts of the flare region* 

Kane, Pick and Rioult (1980), in atteiipting to explain their inferred 
spatial separation of the individual hard X*ray bursts, rule wt chance coin- 
cidence of the tvfo bursts occurring independently as being statistically im- 
probable^, and point out that any signal from one region to the other would 


^ Kane, Pick, and Racult (1980) quote an occurrence rate of * 10 day**^ for the 
class of burst appropriate to their observations, and so dedi^ce the probabil- 
ity of two bursts occurring within a five second time interval (the temporal 
separation of the bursts) to be < 10 This is in fact the probability, 
based on Poisson statistics, of two bursts occurring i*** a given 5 second 
time interval. However, since the time of the first bui^t is arbitrary, one 
should actually calculate the probability of a single burst occurring in a 
given 5 second interval; this evaluates, using the same occurrence rate 
quoted by Kane, Pick, and Raoult, to be • 5 x 10“^, or about 1 in 2000* 
Althou^ this is still a small number, it is conceivable, considering the 
large number of bursts observed, that chance coincidence may in fact have 
been the responsible agent in the event st\uiied by Kane, Pick, and Raoult. 

have to travel much faster than any reasonable Alfven speed, thus ruling out a 
sympathetic trigger. Ttiey thus interpret their observations as in^lying acti- 
vation of two different source regions by the same external disturbance, i.e. 
an energetic electron beam in a "thin target" (Datlowe and Lin 1973) scenario* 
However, in the event studied by Kane, Pick, and Raoult (1980), the two flar- 
ing regions are separated by a distance confer able to a solar radius, and 
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therefore clcirly belong to difforent activo ragicnai for two hard X^ray 
sources located in tha same active region, a sympathetic trigger cannot be 
ruled out by such disturbance velocity considerations (see Karpen, Crannell, 
and Frost 1979). In this paper we propose a means of producing just such a 
trigger* 

To date, all models of the response of the solar atmosphere to a flare 
energy input have been one«-dimensional in nature* This simplifying character- 
istic is usually justified by noting that for typical parameters in the flar- 
ing corona and chromosphere the quantity 6, defined as the ratio of gas to 
magnetic pressures, is much less than unity, so that the plasma is constrained 
to move along the magnetic field lines* (The a88ung>tion ofB < 1 breaks down 
near the photosphere [number density n > 10^^ cm*^l ,but such deep layers of 
the atmosphere remain relatively undisturbed during the flare [Machado, 

Ensile, and Brown 1978].) However, the relation 0 € 1 follows from considera- 
tion of steady-state conditions; as we shall show below (§11) the fact that 
the timescales for impulsive heating, and for hydrodynamic relaxation, of the 
solaur chromosphere are markedly different can result in the formation of re- 
gions of transiently very large gas pressure, in which B is confer able with 
unity* While to explore fully the resulting three-dimensional time-dependent 
response of the solar atmosphere to the flare energy input is beyond the scope 
and purpose of the present paper, we shall discuss scnni -quantitatively the ef- 
fect of creating such a transient high- 6 regime. In particular, we shall dem- 
onstrate (§111) how this high gas pressure region causes the energized flare 
loop to expand laterally, dragging the (frozen-in) magnetic field lines with 
the plasma* A magnetohydrodynamic disturbance is thus established, which can 
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subsequMtly interact with thm fittid lines of a nsighborinq quasi^static loop 
in the same solar active region# causing magnetic reconnection and further re- 
lease of energy# thus triggering the next EFB* 'niis process can be continued 
amongst the various loops in the preflare active region cooplex# producing the 
"multiple-spike” X-ray flux versus time profile characteristic of so many 
events (Hoyng# Bro%m# ar1 van Beek 1976; Karpen, Crannell# and Frost 1979). 

The energetics and timescales of this process are explored in fill# and it is 
shown that both are coi^atible with their counterparts inferred from hard X- 
ray data* In fiv we state our conclusions and suggest areis for more detailed 
research* 


ll* FORMATION OF REGIONS OF TRANSIENTLY HIGH GAS PRESSURE 
DURING THE IMPULSIVE HEATItK? PHASE OF A FLARE BURST 


There is some controversy as to the exact mechanism responsible for heat- 
ing ths solar chromosphere during flares# candidates including non-thermal 

electron bombardment (Bro%m 1973; Lin and Hudson 1976; Emslie 1980)# soft X- 

* 

ray irradiation (Somov 1975; Henoux and Nakagawa 1977, 1970; Machado 1978; 
Hyder 1981)# and dissipation of hydrodynaunic (Craig and McClymont 1976) and 
thermal (Smith and Harmony 1981) shocks* For the purposes of the present dis- 
cussion we shall concentrate on the first of these# for two reasons. First# 
our discussion will relate to events with a strong hard X-ray signature; to 
date all models (both "thermal" and "non-thermal"; see Emslie and Rast 1979) 
of the hard X-ray production mechanism involve a substantial amount of the 
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flat* anarv to b« in tha font of a baaa of pracipitatin? high anargy alac* 
trons (saa, a»g«» Brown 1971 i Malroaa and Brown 1976i Brown« Nalroaa# and 
Spicar 1979) Vlahoa and Papadopouloa 1979) Cmalia and Vlahoa 1980) Smith and 
Brown 1980). Sacond, and mora importantly for our praaant study, this macha* 
nimn of ehromospharic hasting has baan axtansivaly studiad in tha litaratura) 
in particular, calculations of tha datailad tima-dapandant rasponaa of tha 
solar atnoaphsra to such a baam of non>tharmal alactrons hava baan carriad out 
by a numbar of authors (Kostyuk and Pikal'nar 1975) Sarnulina at al . 1980) 
Somov, Syrovatskll, and Spaktor 1981). 

Early work on the response of the solar chromosphere to a hypothetical 
flare energy ii.put %fas, for reasons of simplicity and tractability, confined 
to a discussion of steady -state model atmospheres (eag#^ Brown 1973; Henoux 
and Nakagawa 1977, 1978; Brown Canfield, and Robertson 1978; LaBonte 1978) a 
Xostyuk and Pikel'ner (1975) attempted a time-dependent calculation of the 
hydrodynamic response of the solar atmosphere to a sustained (« 100 s) injec- 
tion of non-thermal electrons a However, they did not adequately treat the 
problem of radiative instability of the heated plasma; further, their pub- 
lished results lack important information, such as density profiles (see 
Canfield et al a 1980) a Recently more detailed work on the time-dependent re- 
sponse of the atmosphere to an injected beam of non-thermal electrons has been 
carried out (Sermulina et ala 1980; Somov, Syrovatskii, and Spektor 1981). 

Due to numerical difficulties, these analyses are restricted to only short 
< 10 second) bursts of heating, and subsequent relaxation of the energized at- 
mosphere to its pre-flare state* This timescale is, however, quite acceptable 
for studying events composed of discrete EFB*s (de Jager and de Jonge 1978; 
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K«rp«n, Crann«ll, and Frost 1979). Furthar. tha rasults claarly show that in 
all casas tha alactron tanq^ratura of tha praflara chrosnspharic natarial Is 
haatad to coronal tanparaturas in tiaws on 'tha ordar of ona saeond. whila tha 
density of this natarial stays ralativaly constant, dua to tha such longer hy» 
drodynamic tine scale. Knowledge of tha detailed behavior of tha atnosphara 
long after this isg>ulsive rise in electron tea^rature is not necessary for 
the subsequent discussion. 

Figure 1 Figure 1 (fron Somov, Syrovatskii, and Spektor 1981) shows the density 

n (cm**^) and electron teng>erature (K) as a function of colunn density N 
-2 

(cm ) after 5 seconds of heating by a beam of non^thermal electrons with an 

—2 —1 

injected energy flux spectrum (electrons cm s per unit energy) in the form 
of a truncated power law: 

V'o> - 21 . »„>>,. m 

®1 ’ 

where 6 ■ 3, ■ 10 keV, and 5 (the total injected energy flux) ■ 1o'^ ergs 

cm 8 • The behaviors of n and in the (hydrostatic) preflare atmosphere 

are shown dasheda Hote that in this preflare state is constant ( ■ 6700 K) 
and n, consequently# has an exponential variation with height (so that n « H). 
This initial state is admittedly somewhat unphysical; however# it turns out 
(Somov# Syrovatskii# and Spektor 1981) that the detailed structure of the 
preflare atmosphere at low column densitites (great heights) is not an impor- 
tant consideration# This is because at such low column densities the flare 
energy input is much larger than any preflare energy source or sink (see 
Emslie 1980)# so that the plasma is explosively heated# resulting in a final 
structure which is essentially independent of the assumed initial 
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conditions* At grsstsr coluim dsnsitias (balow tha ragion of axplosiva 
haatin 9 — taa Brown 1973), tha praflara and flora anar^y inputs ara nora 
ccoparabla, and dua considaration mist ba givan to tha initial conditions of 
tha atnosphara, Sonov, Syrovatskii, and Sjpaktor's (1961) assumad valua of 
- 6700 K raflacts this point, sinca this is a typical taiqiaratura for tha 
uppar chronosphara (saa Machado at al . i960). Claarly a mora datailad 
dascription of T^(N) for larga N %rould ba dasirabla in ordar to mora accurata-’ 
ly assass tha structura of tha flaring chromosphara at thasa lavala, Howavar, 
this papar will ba concarnad only with tha structura of tha atnosphara in tha 

axplosivaly haatad ragion, and so wa shall not dwall on this mattar furthar. 

Also sho%m in Figura 1 is tha flaring prassura 

P • nT^ cn”^ K (2) 

and tha equipartition magnetic field strength 

• (8irkP)'/2 G , (3) 

where k is Boltzmann’s constant. In writing eq. (2)/ %<e have neglected the 

contribution from the ions, since (Somov, Syrovatskii, and ^ektor 

1981). It is clear from the figure that ^ region of greatly enhanced gas 

pressure is formed in the range 10^^ < N < 4 x 10^^ cm"^. For 10^^ 1 ^ $ 2 x 

10^^ cm"^ this is the result of a sudden increase in temperature over pre flare 

19 

values (the density remaining essentially unchanged), while for 2 x 10 
cra”^ < N < 4 X 10^® cm”^ tha excess pressure is attributable to a density en- 
hancement, brought upon by compression due to the overlying high pressure 
material. This ragion of enhanced density is vary thin, being only dN/n > 2 


kn thick, and is a characteristic feature of the chromospheric side of a 


(transient) flar«-»induc«d transition rsgien* Whsn hsatsd it sxpandt at con- 
stant prsssurs to form a region with the original prsflars density and a new 
hotter, temperature, and a new region of enhanced density is created so mewh at 
deeper in the atmosphere* As the heating proceeds, this layer (and so the 
transition region) "burns** its way deeper and deeper into the atmosphere* 

This process stops when the heat being deposited on the chromospheric side of 
this region has been reduced, by attenuation of the electron flux by the 
overlying material, to such an extent that it can no longer overcome the 
ability of the plasma to radiate the energy away* The plasma is now unable to 
enter the temperature regime of radiative instability (see Field 1965; >x and 
Tucker 1969) and so heat up to coronal teiig>eratures (* 10^ K)* The atmosphere 
thus remains approximately in this state (ignoring heat redistribution due to 
thermal conduction and radiative transfer) until the cessation of the electron 
energy input, at which point it starts to relax back, through both thermal 
conduction and radiation, to its preflare state. 

Although the results of Somov, Syrovatskii, and S|)ektor (1981) are rather 
specific, ym may use the physical insist gained in the above discussion to 
generalize their results to other cases. The level of the transition 
region will be determined by the total amount of energy available for heating 
material to coronal temperatures, while the magnitude of the gas pressure in 
the transient regime at this newly-formed transition region will be approxi- 


mately given by 
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wh«r« 9^ is th« initial prstturs at M - Tg And coronal and 

chrMoas^aric alactron tafnparaturaa (• lO^K an< 10^ K) raapactivaly# and tha 
factor 2 in tha danoodnator follow from tha fact that during tha im* 

pulaiva phaaa of tha flara^ whila in tha praflara chronoapharas 

Wa now procaad to aatimata from ai^la thaoratical conaidarationa* 

Tha haating rata (arga ca*^a*^) due to Coulooib colliaiona of a baam of non- 
thamal alactrona, injactad vartically with anargy apactrum T^(Z^) ia (Bro%m 
1973; toalia 1978*-w ignora tha affact of ravaraa currant ohmic loaaaa in tha 
impulaiva phaaa [aaa &nalia 1980 ^ 1981] ) 


Ig(N) - Kn / 


F (E )dE 
o o o 


E* E (1-31W/E^)^'^^ 

o o 


(5) 


4 

whara K ■ 2ira A (a « alactronic charga [a.asu*]# A • Coulomb logarithm), and 


B* • max (e^, [ 301 ]’/^) , 


( 6 ) 


baing tha low anargy cutoff to tha injactad baam wnargy apactrum— aaa aq. 

—2 — 1 

(1)« The total heating rata (ergs cm ^s ) do«m to column depth N is thus 

N H I (N')dN' 

J,(») - / I («•) ji. d»' - / ;; , (7) 

0 0 

or, uaing aqa. (1) and (5), 


- K (6-2) 5 


6-2 


N 


/ / 


dE dN* 
o 


N'-O E -E* E *”^’(1-3KNVE^ 
o o o 


( 8 ) 


For N < > l^/3K, m h«v« I* ■ ravarfing th« ordar of intagration than 

/ialda 


J,(H)| • f [l 




(1-3KS/*^)''^*dl ] . 
o o 


(9) 


Sattlng aac raducaa to a foni convanlant for computation: 


Jg(H) 


I - f [l - 

N<N. 


(«- 2 ) 


J CO* 0 


(1 - 


2 flvV 3 , 

V coa 0) sin 


0 d0] 


( 10 ) 


where 


V 


3W 

‘2 


E 



( 11 ) 


For N E* ■ (31Cm/‘^ 2; thue, in evaluating for N>N^ we must split the 

integral over N’ into t%a> parts— one from N*»0toN*»N^ and the other from 
N* • to N* • Na The first part, is given by equation (10) with v • 1# 

viz. 

j b(| , 3)] , (12) 


wh.ra B is tha Bata function (saa Ahranowlts and Stagun 1965, p. 258). the 
second tarn. 


Jb2<N) 


(«-2)5E 


6-2 


dE dN* 
o 


’ N'-N E ■(3IW')''^2 (i-3»WVE^ 

I o o o 




(13) 
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can be evaluated in a straightforward manner to give (see Bmsllef Brown, and 
Donnelly 1978) 


“I® y 


(14) 


Finally, adding equations (12) and (14) gives 


j (i'j)i '■ 

® N>N^ 3 2 3 1 


-f 


(15) 


note that this result can also be obtained by using the identity 


« I_(N' )dN’ 
} - J^ci . J3(N) ♦ 


(16) 


where the second term is evaluated in a strai^tforward manner with E* = 
( 3KN • )^2 . 


Equations (10) and (15) are the results needed to continue our analysise 
In the case 6 = 4 (a not unreasonable value— -see Hoyng, Brovm, and van Beek 
1976) they reduce to the simple expressions 


JjIN) 


3N 

y(i-— ) ' 


We shall use these expressions hereafter. The effect of varying 


(17) 

6 will be 




considered briefly belcr^ 
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Now consider & beam injected for a tine t, which deposits its energy 
above column depth Hp|^ and heats this material from chromospheric to coronal 
temperatures* Energy balance in the region M < thus dictates that 


T - N™„h(T -T . ) “ N kT , 
B TR TR c ch TR c 


(18) 


where we have Ignored energy losses due to thermal conduction and radiation (cf. above) 
this neglect will cause the calculated from eg. (18) to be an upper limit. 

IQ .2 

For » 10 keV (see Somov, Syrovatskii, and Spektor 1981), - 10 cm 

(Eraslie 1978); thus, by appealing to Somov, £/rovatskii, and %>ektor*s (1981) 
results, we may take N,pj^ ^ ^1 shall in fact verify this assumption ^ 
posteri ori ) e Using eq# (17) in eq# (18) then yields 


1 




(19) 


where 


4R 


3N 


N.kT 

TR ^ 1 C 

; C 


1 


JT 


Equation (19) solves to give 




X » [l + (1-3^)^2 ] ^ 


( 20 ) 


( 21 ) 


where the positive sign in the square brackets has been chosen, so that as 

(^-^O), > N.e In terms of physical variables, eq. (21) is 

TR 1 
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“t»‘ wT I’ * • 

c 


( 22 ) 


Now, by eq« (4), the transient (i.e* flaring) pressure at this level is given 


by 


N„_m„g T 
TR H^e c 


2)cT 


ch 


(23) 


where mjj is the hydrogen mass and the solar gravity. Substituting for 
from eq. (22) and inserting numerical values gives 

« 41.4E^T i, 

P » 6 X 10^’ (l+[l ^2 } . (24) 

ch 


We may now compare this expression with the results of Somov, 
wyrovatsKii, and Spektor (1981); see Figure 1. Ignoring the discrepancy in 
the value of 6 used (for other 6 in the range 3-6 the integrals in eqs. [10] 
and [13] may be evaluated numerically [or analytically] to yield results which 
differ from those for 6 » 4 by only some 10%), we find from eq. (22), with J « 
10^^ ergs cm"^s"\ * 10 keV, x * 5s, =* 10^ K, and * 10^ K, that 


TR 


» 1.8 X 


4«20 -2 
10 cm 


(25) 


(note that this is indeed greater than N^, as was earlier assumed). Further, 


eq. (24) gives 


19 -3 

P =» 6 X 10 cm K 


( 26 ) 
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corresponding to (eq. [3] ) 


B « 450 G . 

eq 


( 27 ) 


Although these values are somewhat larger than the values given by Somov# 
Syrovatskii# v^nd Spektor (1981) (see Figure 1), we note that the value of 
in eq« (25) is more coR^tible with empirical estimates of this quantity 
(Machado et al> 1980), possibly indicating deficiencies in the calculations of 
Somov# Syrovatskii# and Spektor (Emslie# Brown# and Machado 1981)* Even allow- 
ing for the fact that our estimates of are necessarily upper limits# due 
to the oversimplistic energy balance equation (18)# it nevertheless appears# 
therefore, that the gas pressure P can easily attain transient values of order 
10 cm K during the impulsive phase of flares# correponding to an equipar- 

tition magnetic field strength 200 G. (For example# increasing T to 10s 

“4 

in the calculations of Somov# Syrovatskii# and ^ektor [1981] should increase 

Vo 

^ ^ factor " 2 ^ over the values in Figiire 1.) 

eq 

These values of B^^ are not small compared to the ambient field strength 
in the chromosphere* Thus it appears that the assuirg^tion of a plasma con- 
strained to move along preexisting field lines is not adequate in describing 
the impulsive phase of a solar flare burst# since there is no longer an over- 
whelming force (due to magnetic pressure) restricting motion perpendicular to 
these field lines* In the next section we shai . discuss the implications of 
this situation for the evolution of the flare. 
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igure 


III. THE INTERACTK^ OT NEIGHBORING LOOPS 


Consider now the situation depicted in the left sketch of Figure 2# which 
shows two typical loops in an active region loop complex with a simple bipolar 
magnetic field topology* (For clarity# the field lines between and around 
loops L^ and L 2 are not shown*) At scxcie instant loop L^ becomes unstable (due 
to# for exaunple# the formation of multiply connected magnetic islands in the 
loop [Spicer 1976# 1977] or the emergence of new flux at its base [Heyvaerts# 
Priest# and Aist 1977] ) and flares# as indicated by the star at the loop apex* 
This causes a burst of hard X«ray emission to occur* Energy is transported 
along the loop and deposited at the dense footpoints# causing a transient 
strong pressure buildup ($11)* Since the solar atmosphere is highly conduct* 
ing# the field and plasma are "frozen together# and so the large pressure 
gradient between the inside and outside of drives a lateral expansion of 
the field lines defining loop L^. (Of course# there are in addition motionti 
along loop driven by longitudinal pressure gradients [see Craig and 
McClymont 1976].) Although a full three-dimensional magnetohydrodynamic 
treatment of this expansion process is beyond the scope of the present paper 
(see §IV for further discussion)# to order of magnitude one may assume that a 
disturbance of the surrounding field lines proceeds horizontally at the ion 
sound speed c^ * * which is approximately equal to the Alfven speed 

« B/(4irnm^)^2 since the gas and magnetic pressures are comparable (Fig* 
ure 1 ) • After a time t the disturbance encounters the neighboring loop L 2 
(Figure 2)# allowing L 2 to release its stored free ener^ (in the form of 
currents) through interaction with the disturbed field lines in a 
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n«utral-she«t>typ* configuration (sea Sturrock 1968), and giving rise to 
another burst of hard X*ray emission, spatially distinct from the location of 
the first burst (see fl). 

In order for the phenomenon just described to explain the features of 
EPB's, the interaction time r must be roughly equal to the temporal separation 
^EFB EPB's* Thus, the lateral separation D of the loops and I <2 must be 
of order 


'^A^EFB 


(28) 


idilch, setting ■ 5 - 10 s (de Jager and de Jonge 1978) and Vj^ " 5 x 10^ 

cm s~^ (using B ■■ 500 G and the density values shown in Figure 1), gives 


D - (3 


5) X 10® 


cm , 


(29) 


an entirely reasonable value* 

The situation envisaged (” driven merging flux'*) is somewhat similar to 
that in the flare model by Gold and Hoyle (I960), or in the emerging flux 
model of Keyvaerts, Priest, and l^st (1977), However, there are some signif- 
icant differences, as Figure 2 shows* Because the two loops under considera- 
tion are part of the same active region complex, which we assume to have a 
simple bipolar structure, the polarities of the corresponding footpoints of 
the two loops are the same ; such a situation gives no reversal of the toroidal 
(i.e. longitudinal) component of the magnetic field at the contact plane be- 
tween and L 2 and so this toroidal component of the field cannot be respon- 
sible for any subsequent energy release as the model of Heyvaerts, Priest, 
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and Ruat (1977). Purthar, tha (pololdal) currants corraapondlng to thaaa tor- 
oidal flalda ara antlparallal in tha contact plana and tharafora rapal , in- 
staad of attracting in tha (fiald-ravaraad) Gold and Hoyla (1960) modal. How- 
ever, by tha same conaidarationa, any twiat (i.a., departura from a potantial 
configuration) of tha magnatlc field linea will be in tha same sanaa as viewed 
by a remote observer, which leads to a reversal of the pololdal cong>onent8 of 
the magnetic field in the contact plane, providing a situation suitable for 
reconnection of this pololdal congxjnent in a neutral-sheet-type scenario (see 
Sturrock 1968). The force required to overcome the repulsion of the pololdal 
currents is the hydrodynamic pressure of the heated transition region in the 
flaring loop; a simple calculation sho«ra that the condition 6 ■ 1 is equiva- 
lent to the condition that the gas pressure gradient can overcome the electri- 
cal repulsion of these pololdal currents. 

This "pololdal flux annihilation" aspect of the present model gives, upon 
further consideration, some rather pleasing results. First, since the toroi- 
dal component of the magnetic field remains intact after the reconnection pro- 
cess, the basic spatial configuration of the loops remain intact, permitting 
the formation of the soft X-ray emitting post-flare loops frequently observed 
(e.g., Pallavicini, Serio, and Vaiana 1977). Second, during the reconnection 
3B/3t, and so V x E, is in the pololdal direction. This gives an induced E 
vector which is predaninantly toroidal, i.e., parallel to the remaining mag- 
netic field lines, thus allowing efficient acceleration of particles along 

these field lines. (In standard neutral-sheet-type models (e.g., Sturrock 

1968] 3B/3t and Y 1 ™®inly vertical, resulting in an induced E vector 

which is horizontal, i.e., perpendicular to the remaining [closed] field 


19 


Iin«8«) Third, there is no need for a complex magnetic topology at the loop 

footpoints, such as is required by both the Gold and Hoyle (I960) model and 

the emerging flux model of Neyvaerts, Priest, and Itist (1977) • 

One may also estimate the energy released in the annihilation of this 

poloidal ^lux* Assuming a twist of » w in each flux tube (see Barnes and 

Sturrock 1972), a toroidal field strength of some 300 G (see above), and a 

2 17 2 

footpoint area A*D *3x10 cm ,we find that the total energy released 
is 


2 2 2 

€ B D L fD\2 
®rel * Sir ' 


(30) 


where L is the length of the flux tube (• 3 x 10^ cm) and e is the fraction of 
the poloidal field that is destroyed. Now.* the proposed model relies on the 
fact that e < 1, since some poloidal flux must r«nain after each loop **encoun- 
ter" if the "domino" effect leading to successive EFB*s is to work. Thus, 

Q Q 

setting e ■ 1/3, B - 500 G, D • 5 x 10 cm, and L ■ 3 x 10 cm, we obtain 


"rel 


» 3 X 


10 ergs 


(31) 


We may compare this with the energy required to produce a single hard X-ray 
EFB. Hard X-ray spectra in moderately large flares may be represented well by 
the power-law form 


-Y -2 -1 -1 

1(e) ■» ae photons cm s keV , 


(32) 


with a ■ 10 ^md Y(*6-1) *3-4 (see Hoyng, Brown, and van Beek 1976). To 
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produc* «ueh a burst by thick target bramsstrahlung of a baaa of high energy 
electrons requires an Injection rate of electrons with energies > 25 keV of 


Fjg “ 4 X 10^^a(Y“1)^B(YJ>^ ,V 2 )(25)"^ - electrons s"^ (33) 


(Hoyng« Bro%m, and van Beek 1976), corresponding to a total Injected power 

P - X 25 keV » 3 X 10^^ - 3 x 10^® erg s“^ . (34) 

Comparing eqse (31) and (34) with a FVHM for the EF6 of *58 (de Jager and de 
Jonge 1978) indicates agreement to within an order of magnitude, which consid- 
ering the crudeness of the above argument, is quite acceptable. 

This completes the scenario of the interacting loop model, whose energet- 
ics and timescales have been shown to agree well with observations. In the 
final section vie shall discuss some features of the model in greater detail, 
and also point out areas for future work. 


IVe DISCUSSION 


The analysis of the preceding two sections is of course somewhat schema- 
tic, and is intended only to provide a general description of the main fea- 
tures of the model. Nevertheless, it not only accounts for the EFB descrip- 
tion of hard X-ray bursts, but also, for plausible parameters, successfully 
accounts for both the observed intensity and tes^ral spacing of ETB's. The 
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only • priori rogulrwaonf of th« modol at* (i) th« initial anargy ralaaaa in 
ona of tha activa ragion loops in Flgura 2), which may occur by any of tha 
currantly propoaad aachaniams (a>g>f tha taaring mods instability in a 
strassad arch gaomatry“~Splcar 1976» 1977i van Hovan 1979), and (ii) tha axis** 
tenca of a larga seals currant pattam in tha activa ragion, to provide tha 
pololdal flux cong>onant which ia annihilatad in tha "domino" procaas of auc- 
cesslva energy ralaases. 

Clearly many aspects of the modal require more detailed study* For exam~ 
pie, in Figure 2 it is tacitly assumed that tha energized flux tube expands 
along Its entire length as a result of the gas pressure forces which are 
localized at its base. Our intuitive justification for this is that the ten- 
sion in the siaignotic lines of force will tend to "straighten out" any "bulge" 
formed by a local expansion of the field; however, the timescale for the en- 
tire loop to respond to this tensile force will be of order corona' *'*’^®*' 

for very large (L/D) may be larger than the inter-loop collision timescale 
(eq. [29]), even allowing for the reduced density, and corresponding higher 
Alfven velocity, in the corona. A full three-dimensional magnetohydrodynamic 
analysis of the situation is clearly necessary to resolve this issue satisfac- 
torily, just as it is also required to adequately model the expansion of the 
flux tube in the (difficult) 6 • 1 regime (see remarks at the beginning of 

In addition to these issues, the details of the Interaction between "col- 
liding" flux tubes (see Figure 2) needs to be explored more fully. This re- 
quires study of driven magnetic reconnection processes, such as have been dis- 
cussed previously, for the case of flux driven by magnetic buoyancy , by 
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H«yv««rt8# Pri«st^ and Raat (1977), and by rafarancaa tharain. This will an* 
abla ua to battar datarmina tha taaporal and spatial charactaristics of tha 
anargy ralaasa in tha quiascant loop (loop L 2 in Figura 2). 

I thank J. Lsach for stiimilating discussions which lad to tha idaas da- 
valopad in this papar. A* N. HcClymont for usaful discussions ragarding tha 
hydrodynamics of tha solar atmosf^ara during flaras, and S* K« Antiochos and 
P« A. Sturrock for thair commants on tha manuscript* Tha author has banafitad 
from discussions hald at tha Study of fiiargy Relaasa in Flarac (SERF) Workshop 
held at Stanford in August# 1980i this %forkshop %#as sponsored by tha National 
Science Foundation and by tha Scientific Committaa on Solar-Tarrastrial 
Physics (SCOSTEP)* Financial support from grants NASA NGL 05-020-272# NASA 
NAGW-92# and contract ONR NOOO 14-75-^-0673, is also gratefully acknowledged* 
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FIGURE CAPTIONS 


Fig* 1* - Variation of density n (cm“^) and electron temperature (K) with 
column density N (cm )# from the Somov, Syrovatskii, and %>ektor (1981) 
calculations of the hydrodynamic response of the solar chromosphere to an 
energy input in the form of a beam of non-therraal electrons* The preflare 
values of n and T^ are shown dashed; see §II for a discussion of these initial 
conditions* The values in the flare correspond to a time 5 seconds after the 
start of the heating and correspond to the largest enhancements in their 
calculations* The ion temperature (not shown) is much less than T^. Also 
shown is the pressure P » nT^ (cm*^ K) and the equipartition magnetic field 
®eq ®®® [3]). Note the large values of P and around the 

transition zone, caused by the in^ulsive heating of plasma; the value of B 

eq 

there is of order 100 G* 


Fig. 2. - Schematic of the flux tube interaction. Loops and L2 are typical 
members of the same bipolau: active region complex and lie across the neutral 
line (heavy dashed) as shown. Their magnetic fields are also twisted in the 
same sense due to a global active region current field J. At some time t * 0 
loop flares (the energy release being signified by the star at the loop 
apex), causing energy to flow downwards toward the chromosphere (the energy 
transport mechanism is here taken to be aon-thermal electron bombardment, so 
that we may analyze its effects in some detail - see §II - however, other 
transport mechanisms [see §11] will have the same qualitative effect on the 
evolution of the system). Upon interaction with the chromosphere, the trans- 
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ported energy heats the plasma to very high (« 10^ K) tengmratures, causing a 
large pressure enhancement at this point (see Figure 1). This causes loop 
to expand# driving a horizontal nagnetohydrodynamic disturbance (small hori- 
zontal arrows at rhe loop base)* At some later time the surface of the 
disturbed field lines has the appearance (right figure; the original loop 
is sho%m dashed)# and is in contact with loop L 2 # with the poloidal com- 

I 

ponents of the magnetic fields of and 1*2 anti.parallel in the contact 
plane* These components reconnect in a nautral-sheet-type scenario and cause 
an impulsive energy release in the hitherto quiescent loop 

For the purpose of the Illustration# the poloidal component of the mag- 
netic fields has been exaggerated* Also# loop has been constructed as a 
simple radial dilation of loop this is not necessarily justified in prac- 

tice (see §IV)* 
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